DETECTION OF C3 IN DIFFUSE INTERSTELLAR CLOUDS 
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ABSTRACT 

The smallest polyatomic carbon chain, C3, has been identified in interstellar clouds 
(At, ~1 mag) towards C Ophiuchi, 20 Aquilae, and C, Persei by detection of the origin 
band in its ^^n^ — electronic transition, near 4052A. Individual rotational lines 

were resolved up to J=30 enabling the rotational level column densities and temperature 
distributions to be determined. The inferred limits for the total column densities (~1 
to 2x10^^ cm~^) offer a strong incentive to laboratory and astrophysical searches for 
the longer carbon chains. Concurrent searches for C^, C2 ^"^^ ^^re negative but 
provide sensitive estimates for their maximum column densities. 

Subject headings: ISM: molecules-Cs , C3 , C^, C2 



1. Introduction 

Carbon chains play a central role in the chemistry and spectroscopy of interstellar space. The 
detection of cyanopolyacetylenes in dense interstellar clouds (Kroto 1981) led to the suggestion by 



^Visiting Astronomer, Canada-France-Hawaii Telescope, operated by the National Research Council of Canada, 
the Centre National de la Recherche Scientifique of France, and the University of Hawaii. 



- 2 - 



Douglas (1977) that carbon chain species be considered as candidates for the diffuse interstellar 
bands (DIBs), which are found in the 4000-8500A spectral region of stars reddened by interstellar 
dust (Herbig 1995). Since then, many other molecules with a carbon chain backbone have been 
identified at radio frequencies in dense clouds (Kawaguchi et al. 1995). Meanwhile, advances in 
laboratory measurements have provided an understanding of the types and sizes of carbon chains 
which have strong electronic transitions in the DIB range (Maier 1998). It is thus somewhat 
surprising that as yet among the bare carbon species only diatomic C2 has been identified in 
interstellar clouds where DIB are detected. 

The 4052A electronic band system of C3 was first detected in comets (Swings 1942) and then in 
circumstellar shells by infrared spectroscopy (Hinkle et al. 1988). Most recently C3 was identified in 
a dense cloud using sub-mm measurements of its low frequency bending mode and Ntot(C3)'^ 10^^ 
cm~^ was estimated (Cernicharo et al. 2000). Snow et al. (1988) established an upper limit of 5x 10^° 
cm~^ for the column density of C3 in the direction of C Oph, some two orders of magnitude lower 
than that set by Clegg k. Lambert (1982). Haffner &; Meyer (1995) made a tentative detection of 
C3 towards an eighth magnitude star in the same part of the sky, HD 147889, at a column density 
of 4x10^^ cm~^. Unfortunately, the star turned out to be a (previously unknown) double-lined 
spectroscopic binary which limited their sensitivity. 

This letter presents the detection of C3 towards three stars and infers the column densities 
in the diffuse clouds. Although C^, and C3 were not detected, upper limits for their column 
densities are estimated. Diatomic species, such as CH, CN, C2, and CH"'", have been detected 
towards two of the stars chosen and their column densities are considered standards with which to 
compare models for the physical and chemical processes in diffuse regions (van Dishoeck & Black 
1986). The four bare carbon chains, C3, Cg", C2 and C3 , were selected for the present search 
because gas phase electronic transitions for these species have been identified in the laboratory in 
the 4000-5500A region and their oscillator strengths are known (see table 3). 



2. The Observations 

Observations of the reddened stars C Oph (HD 149757), C Per (HD 24398) and 20 Aql (HD 
179406) were made with the Gecko echellette spectrograph on 2000 July 16 and 19, fiber fed from 
the Cassegrain focus of the Canada-Prance-Hawaii 3.6-m telescope (CFHT) (Baudrand k. Vitry 
2000). All three stars have a visual extinction, A^, near 1 and were chosen because they are bright 
with sharp interstellar K I lines indicating either single clouds or little Doppler distortion (in the 
case of C, Oph, ? resolved the C2 at 8756 A into two close velocity components separated by 1.1 
km s"-*^). 

The detector was a rear illuminated EEVl CCD (13.5 jxm^ pixels) and the spectral regions 
were centered at 4047 A in the 14th order, and at 5060 A and 5400 A in the 11th and 10th orders, 
respectively. The ultraviolet Gecko prism was used to isolate the 14th order, the blue grism for the 
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11th order, while the stock CFHT filter #1515 was used for the 10th order observations. Individual 

spectra had exposure times ranging from 5 to 20 minutes and were obtained with continuous fiber 
agitation to overcome modal noise. The resulting combined spectra for the individual stars at each 
wavelength had unusually high signal-to-noise ratios (S/N~800-4000) for CCD observations. The 
Th/Ar comparison arc spectra, taken before and after each spectrograph wavelength reconfigura- 
tion, had a typical FWHM of 2.8 pixels, which corresponds to resolutions of i? = 121000, 113000 
and 101000 at 4047, 5060 and 5400 A, respectively. 

Processing of the spectra was conventional. Groups of biases were taken several times through- 
out each night and at each grating setting a series of flat-field spectra of a quartz-iodide lamp were 
recorded. The biases and appropriate fiats were averaged and used to remove the zero-level offset 
and pixel-to-pixel sensitivity variations of the CCD. One-dimensional spectra were extracted us- 
ing standard IRAF routines. Spectra of Vega and other unreddened stars were used to search for 
contaminating telluric water vapor lines and stellar photospheric features. Heliocentric corrections 
were applied to each spectrum. 

The observations are summarised in Table 1 which lists exposure times and S/N per pixel for 
each spectral region. The final column gives the radial velocities measured from the interstellar 
K I 4044.1 and 4047.2 A lines. These velocities have been applied to each spectrum to put the 
interstellar features on a laboratory scale. 

3. The Results 

The ^^n„ — X^T,'^ origin band of C3 is quite clearly detected towards all three stars. Figure 1 
compares the observed spectra with a simulated C3 absorption spectrum based on the spectrograph 
resolution and assuming a Boltzmann distribution in the ground state rotational levels with a tem- 
perature of 80 K. The continuum noise level in the observations is ~0.1 %. Low order polynomials 
have been applied to the stellar data to give a level continuum (base line) and, in the case of C, 
Per, a weak, broad stellar feature at 4053. 2A has been removed. Residual broad features in the 
final spectra are only a few tenths of a percent deep, much less than in the original, and they in no 
way mask the sharp C3 lines. In the simulation the rotational line intensities were calculated using 
the Honl-London factors, while the line positions were taken from the laboratory measurements 
(this avoids the problem of a perturbation affecting the low J ground state levels which is not 
accounted for by the fitted spectroscopic constants) . The individual rotational P, Q and R lines are 
clearly resolved in the spectra of all three stars. Table 2 lists the observed positions and equivalent 
widths of each rotational line assigned in the spectra of the three stars. The table also gives the 
corresponding positions measured for these transitions in the laboratory (Gausset et al. 1965). The 
positions of some 30 lines in the spectrum of ( Oph agree with the laboratory data to within 0.1 
cm~^ providing an unambiguous identification of C3 in these diffuse clouds. 

Figures 2, 3, and 4 show the results of equivalent searches for origin band of C^ at 
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5078.1A (C Oph only), the B^T,+ -X'^Y.+ origin band of at 5415.9A and the A^A^-X'^Ug origin 
band of at 4040. 4A, together with simulated spectra for these transitions at 80 K based upon 
the published spectroscopic constants (Maier & Rosslein (1988), Herzberg & Lagerqvist (1968), 
Tulej et al. (2000)). In the cases of C2 and the linewidth was assumed to be determined by 
the spectrograph resolution. For C3 the excited state has been identified as a short-lived Feshbach 
resonance and the measured natural linewidth of ~1 cm~^ is employed in the simulation. Weak 
telluric lines were been removed from the Cg" and observations using standard procedures. 

Unlike the 4050A region for C3, each of these spectral regions is contaminated by weak stellar 
features. Nonetheless, for the C2" and ions, there are no sharp features corresponding to the 

rotational lines in the simulations. On the other hand, for Cj there are features in the spectrum 
of C Per and (much less convincing) in the magnified plot for Oph which appear to coincide in 
position and shape with the simulated band heads. It is unlikely that these are due to C3 because 
the stellar lines in ( Per have exactly the same shape as the coincident features. For ( Oph the 
whole CJ spectrum sits within a weak stellar feature (the lines are broadened by rapid rotation 
~400 km s~^). The photospheric lines in ^ Oph show nonradial pulsation 'ripples' which will be 
'washed out' to some extent by the long exposure time employed. The spectrum of 20 Aql, which 
normally has the strongest interstellar lines of the three, is free of stellar features but has no features 
coincident with the C3 simulation. It is concluded that, while interstellar C3 might be absent for 
20 Aql and present for the other two stars, it is more likely that in the latter cases the features are 
instead stellar. 

Table 3 gives the measured equivalent widths for the most intense C3 line (Q(8)) for each star 
together with an la error estimate. For Cj, CJ and C3 , 3a detection limits are given. The la 
level errors and detection limits are derived from: 

Wx = iwd)^S/N)-^ 

where the la limiting equivalent width, Wx, and the FWHM of the feature, w, are both measured 
in A, the spectrograph dispersion, d, in A pixel"^, and S/N is the signal to noise per pixel. From 
the simulations, w = 0.045, 1.0, 0.13 and 0.045A for C3, C3 , and C^, respectively. 

In the case of C3, equivalent widths, Wx, were determined for each rotational line (varying 
between 0.3-2. 7x 10~^A) and, in combination with the transition oscillator strength, fo_o, and Honl- 
London factors, the column densities, Nj{Cs), of each rotational level (J) in the ground electronic 
state were calculated (Larsson 1983). In cases where several rotational lines originating from the 
same level were assigned (e.g. P(8), Q(8), R(8)) the mean of the determined column densities was 
taken. Figure 5 shows a Boltzmann plot of ln(A^j/(2J + 1)) vs. the rotational energy (Gausset et 
al. 1965) where the slope is inversely proportional to the rotational temperature. Among the lowest 
rotational levels (J <14) the populations are reasonably approximated by a distribution at 50-70 
K, whereas the higher rotational levels correspond to a temperature of 200-300 K. The simulation 
in Figure 1 uses 80 K as this represents an average temperature for the entire rotational population 
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and allows both the high and low J lines to be identified. The high temperature component of the 
distribution is apparent in the astronomical spectra where the R band head and the higher Q lines 
are more intense than in the simulation (Figure 1). 

Lutz & Crutchcr (1983) also found a bimodal population distribution for C2 in diffuse clouds, 
with similar characteristic temperatures for the low and high J values. The lower temperature 
is interpreted as the kinetic energy of the cloud and for both C2 and C3 the values obtained are 
comparable to those used in models of diffuse clouds (van Dishoeck k, Black 1986). The higher 
temperature component is attributed to repopulation of the levels in the ground electronic state 
by radiative pumping from excited states. In the case of C3 it is expected that both the state 
and the higher lying state will contribute to the radiative pumping. 

The sensitivity of these measurements is such that Nj{Cz) in the range 0.2-2x10-'^^ cm~^ is 
determined for rotational levels up to J=30. The Nj values were summed to give the estimated 
lower limits in the range 1-2x10^^ cm~^ for the total column density, iVtot(C3), in Table 3. 



3.1. Discussion 

A previous search for C3 in the direction of Q Oph did not identify the molecule (Snow et al. 
1988). It is unclear why this was the case as, in the light of the present observations, the signal-to- 
noise quoted for these measurements was adequate and the upper limit given was some thirty times 
lower than the column density reported here. The present measurements for the column densities 
of C3 are of the same order of magnitude as the tentative estimate for a translucent cloud (Haffner 
&; Meyer 1995) and can be compared to those of other polyatomic molecules observed in diffuse 
interstellar clouds. Column densities (also towards C, Oph) in the 10^^""*^^ cm^^ range have been 
inferred for HCO"*", C2II and C3II2 from observations in the mm region by Lucas &l Liszt (2000), 
while Hg^ has been identified in diffuse regions and Ntot(II^) estimated ~ 4x 10^^ cm~^ by McCall et 
al. (1998). Column densities of C2 towards C, Oph and C, Per have been determined in the 2-3x10^^ 
cm~^ range by Federman & Huntress (1989). A current model of the diffuse clouds by Ruffle & 
Herbst (2000) predicts Ntot(C2)/Ntot(C3) ~ 20 (on a 10^ year time scale), implying Ntot(C3)~ 10^^ 
cm~^, in agreement with the values deduced from the Table 3. The main production route to C3 
is presumed to be the dissociative recombination process: CsH"*" + e — > C3 -|- H, where C3H+ is 
produced from smaller species by C+ ion insertion. Under conditions where ultra-violet radiation 
penetrates, photodissociation of C3 takes place at a threshold of 1653A. As the strong ■'^S+— 
electronic transition of C3 is predicted to occur around 1700A (Forney et al. 1996), the dissociation 
process, C3 — > C2 + C, may be an important destruction pathway in diffuse clouds. 

Although only upper limits for the column densities of Cg , Cg" and C^ could be presently 
established, these species are of interest as small ionic carbon fragments play a crucial role in the 
ion molecule schemes for diffuse cloud chemistry (Dalgano & Black 1976). The ion is the only 
bare carbon cation for which the gas phase electronic spectrum is known (Maier & Rosslein 1988). 
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In diffuse clouds it is supposed to be the product of the fundamental step: C"*" + CH — > Cg" + H. 

Its main destruction mechanism is hydrogenation: + H2 ^ C2H+ + H, which dominates over 
recombination with electrons in diffuse regions. The diffuse cloud model (Ruffle & Herbst 2000) 
predicts a abundance a factor of 10~^ lower than C2, implying a column density ~ 10^'' cm~^, 
in accord with the upper limit in Table 3. 

Published models for diffuse regions do not include the smallest pure carbon anions, C2 and 
in their reaction libraries. Unlike C^, C2 does not react with H2 so its main destruction 
mechanism is expected to be photodetachment. The similar rotational line widths and oscillator 
strengths of the and C2" transitions lead to similar upper limits for their total column densities. 
The width of the unresolved bands for C3 and the presence of weak stellar features in this spectral 
region means that a higher column density of this ion could have escaped detection. 

The detection of C3 provides a powerful incentive for the laboratory study of the electronic 
transitions of longer carbon chains in the gas phase with the aim of comparison with DIB data. 
The question as to what types and sizes of carbon chains will have strong transitions in the 4000- 
9OOOA range has already been answered: for example, ^E+— transitions of the odd-number 
bare chains, C2n+i, n=8-30 (Maier 1998). The existence of linear carbon chains up to C21 has 
been confirmed by the observation of their electronic spectra in neon matrices (Wyss et al. 1999). 
As the oscillator strength scales almost linearly with the length of the molecule, one can expect 
fo-o ~ 10-20 for these carbon chains. With such an oscillator strength, a species with a column 
density ~ 10^^ cm~^ would be enough to give rise to a strong DIB, with an equivalent width of lA. 
In view of the column density ~10^^ cm~^ for C3 determined in this work for three diffuse clouds, 
this appears to be a reasonable expectation. 

The support of the Swiss National Science Foundation (project no. 20-055285.98), the Cana- 
dian Natural Sciences and Engineering Research Council and the National Research Council of 
Canada is gratefully acknowledged. The authors thank the staff of the CFHT for their care in 
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Tablo 1. The Observations 



star 


Sp/L 


V 


E(B-V) 
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5400 A 


K I RV 














T S/N 


T 


S/N 


(km s-i) 


C Per 


Bl lb 


2.85 


0.28 


4800 


1200 




2700 


1900 


+13.91 ±0.24 


C Oph 


09.5 V 


2.56 


0.30 


5400 


2400 


5400 4000 


3000 


2200 


-14.53 ±0.18 


20 Aql 


B3 V 


5.36 


0.27 


10800 


800 




8400 


900 


-12.53 ±0.08 



^Exposure times in seconds, 
''per pixel. 
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Table 2. The C3 Rotational Lines 
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59 089 
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3.130 
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52.271 
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52.900 
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53.781 


1.098 
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54.113 
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54.445 
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P(18) 
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''Gausset et al. (1965) 

''These features are assigned to two blended lines of similar intensity in the simulation and are not used to 
determine rotational level column densities. 

R(0) is a very weak line with a correspondingly larger positional error compared to stronger neighboring 
lines. 
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Table 3. Total and rotational level column densities of C3 towards the three stars and estimated 
upper limits for the equivalent widths and column densities of C3 , C2 and 



species 


transition 


Ao-0 


fo-0 


star 












(A) 






lO-^A 


10'° cm-2 


10'2 cm-2 


C3 




4051.5'' 


0.016'= 


COph 


2.3±0.1 


20 


1.6 










20 Aql 


2.2±0.3 


19 


2.0 










C Per 


2.0±0.2 


17 


1.0 






4040.4'^ 


0.04" 


COph 


6.0 




<0.3 










20 Aql 


20 




<1.2 










C Per 


12 




<0.7 


ct 




5066.9* 


0.025' 


COph 


0.35 


<1.1 


<0.04 




S^EJ - X2E+ 


5408.6^ 


0.044'' 


COph 


0.30 


<0.5 


<0.02 










20 Aql 


0.55 


<0.9 


<0.03 










C Per 


0.35 


<0.6 


<0.02 



^For C3 the Nj values refer to the J=8 level and Wx to the most intense Q(8) line. For CJ and 
Cj Wx are 3a detection limits, based on simulated spectra at the spectrograph resolving power, 
and the Nj values are determined from the most intense R lines at 80 K. For Cj individual 
rotational lines are not resolved and Wx is the 3a detection limit for the blended R branch. 

''Gausset et al. (1965) 

■^Becker et al. (1979) 

''Tulej et al. (2000) 

''Estimated value (Tulej et al. 2000) 

'Maier & Rosslein (1988) 

*Herzberg & Lagerqvist (1968) 

''Leutwyler et al. (1982) 
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Wavelength (A) 

Fig. 1. — Rotational lines in the origin band of the C3 ^^H^ — X^T,^ transition towards ( Oph, ^ 
Per and 20 Aql. The bar indicates 0.4% of the continuum. A known, weak stellar absorption line 
at 4053.2 A has b ecu removed from the spectrum of ^ Per. The simulated spectrum is for Trot — 
80 K at a spectral resolution of 110,000. 




I 1 1 1 1 1 

5050 5055 5060 5065 5070 5075 

Wavelength (A) 



Fig. 2. — A spectrum of ( Oph in the region of the C2" -B^S~ — X^E~ origin band. The upper curve 
is a simulated spectrum of C2" showing pairs of rotational lines. The feature at 5062 A cannot be 
stellar as it is too narrow but might be a sharp DIB. Note the very high signal-to-noise (> 4000) 
in the stellar spectrum. 
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5395 5400 5405 5410 



Wavelength (A) 

Fig. 3. — Spectra of C, Oph, ( Per and 20 Aql in the region of the C2 origin -B^S+ — band. 
The upper curve is a simulated spectrum of C2 showing individual rotational lines. The feature 
at 5404.6 A is a known, diffuse interstellar band (DIB). There is a strong, rotationally broadened 
stellar line at 5412 A in ( Oph. 




— I 1 1 1 1 1 1 — 

4034 4036 4038 4040 4042 4044 4046 



Wavelength (A) 

Fig. 4. — Spectra of Oph, ( Per and 20 Aql in the region of the Au — X'^Ug origin band. The 
upper curve is a simulated spectrum of C3 showing two unresolved spin-orbit component bands. 
The spectrum of Per is strongly contaminated by stellar absorption features. The sharp line at 
4044.1 A is the interstellar K I line which was one of two features used to place the spectra on a 
laboratory wavelength. Note that the intensity scale for ( Oph is four times greater than for the 
other two stars. 



25.0 




Fig. 5. — Plot of ln(A^j/(2J+l)) vs. rotational level energy, E(J), of C3 for C, Oph. The error bars 
are ±la in the measured column densities. Straight line fits are shown for J<14 (dashed, 60 K) 
and J>14 (solid, 230 K). 



